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ABSTRACT 



In this paper we discuss the entropy of quantum black holes in the LQG formalism when 
the number of punctures on the horizon is treated as a quantum hair, that is we compute the 
CN ! black hole entropy in the grand canonical (area) ensemble. The entropy is a function of both 
the average area and the average number of punctures and bears little resemblance to the 



C**~ ' Bekenstein-Hawking entropy In the thermodynamic limit, both the "temperature" and the 



chemical potential can be shown to be functions only of the average area per puncture. At 
a fixed temperature, the average number of punctures becomes proportional to the average 
area and we recover the Bekenstein-Hawking area-entropy law to leading order provided that 
the Barbero-Immirzi parameter, 7, is appropriately fixed. This also relates the chemical 
potential to 7. We obtain a sub-leading correction, which differs in signature from that 
obtained in the microcanonical and canonical ensembles in its sign but agrees with earlier 



5-1 . results in the grand canonical ensemble. 
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I. INTRODUCTION 



Loop quantum gravity (LQG) is a background independent non-perturbative canonical 
quantization of general relativity (GR) [Q employing a canonical chart consisting of the 
su(2) Lie algebra valued Ashtekar-Barbero connection |2| and densitized triads as canonical 
variables. Because of the geometric interpretation of these variables it is possible to construct 
operators in the quantum theory corresponding to classical geometric quantities such as 
area and volume. At the kinematic level, these operators have discrete spectra, quantized 
in units of an appropriate power of the Planck length (upto an arbitrary constant, the 
Barbero-Immirizi parameter). At the dynamical level, the question of finding the quantum 
Einstein's equations becomes the problem of finding the state space lying in the kernel of 
the Hamiltonian constraint. Although a solution to this problem remains to be found and 
the full picture of quantum gravity is yet to emerge, the structure of space-time revealed at 
the kinematical level provides an interesting approach to black-hole physics when the black 
hole horizon is modeled as an isolated horizon || . The isolated horizon is taken as an inner 
boundary of space-time, which gets threaded by links of eigenstates of geometric operators, 
namely the "spin- networks" . Each puncture contributes some quantum of area according 
to the su(2) label it carries, which provides a microscopic description of the horizon and its 
degrees of freedom. Standard statistical mechanics techniques may be used to determine 
the entropy associated with the horizon and the results can then be tested against the semi- 
classical theory, in particular to the Bekenstein-Hawking (BH) entropy 0. This program 
has been carried out successfully in LQG. The BH entropy is recovered as a dominant term 
at the macroscopic level and sub-dominant terms of quantum origin are also obtained |5|. [| . 

It is worth noting that the BH area-entropy law has also been obtained in diverse ap- 
proaches to quantum gravity and each of these approaches has succeeded in providing an 
alternative interpretation of the quantum gravitational degrees of freedom comprising a 
black hole. Sub-dominant, logarithmic corrections to the BH law have also been proposed 
in these approaches 0. 

In LQG the black hole is traditionally treated within the microcanonical ensemble, since 
the horizon is taken to be an isolated horizon ||. Use of the microcanonical ensemble is 
generally justified by the fact that because there is no energy flux across the isolated hori- 
zon, the area is expected to remain constant. However, when the possibility of black hole 
evaporation is taken into account, an exchange of the number of punctures or the spin-label 
of a puncture with the bulk geometry can, in principle, give rise to area fluctuations as well 
as fluctuations in the number of punctures. Thus it appears that the canonical and/or grand 
canonical ensembles, in which the number of punctures is treated as a statistical variable, are 
more appropriate choices if one intends eventually to describe physical processes involving 
black hole evaporation. There is another reason why the microcanonical ensemble may not 
be suitable for the study of black hole thermodynamics: it is now well known that the micro- 
canonical entropy, apart from being non-uniquely defined, is not a differentiable function of 
the area but rather a staircase (or ladder) function || |K| . As pointed out in JTT| , this makes 



it difficult to interpret basic thermodynamic variables, such as the temperature, which are 
defined in terms of derivatives of the entropy. The canonical and grand canonical partition 
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functions are assumed, by contrast, to be smooth functions of their arguments. Therefore, 
recently there has been renewed interest in treating the horizon within the canonical and 



grand canonical area ensembles [[T2|, [13| . 

The idea of using the horizon the ensemble variable, as opposed to the more 

traditional energy ensembles, dates back to the work of Krasnov |14|]. It was originally 



motivated by the fact that if a microscopic description of the black hole entropy is to yield 
something like the BH law then the density of states should satisfy a relation of the form 
lnf2(M) ~ M 2 , which would lead to a divergent partition function in the energy ensemble. 



A novel justification for its use may be found in |T3|. The use of the number of punctures 



N, as a statistical variable was advocated by JTbl and recently revived in 12 



In [17] we presented a detailed analysis of black holes with the LQG area spectrum in the 
canonical area ensemble (fixed number of punctures) and obtained a logarithmic correction 
with a negative coefficient, in agreement with the microcanonical entropy results. In [|11|], the 
entropy was computed in the grand canonical ensemble for two counting schemes, viz., the 
one proposed by Domagala and Lewandowski (DL) fl8| and another, more recently proposed 
by Engle, Noui and Perez (ENP) [[B|. The chemical potential was set to zero from the start 
and a sub-dominant logarithmic correction with a positive coefficient was determined in the 
DL case whereas no correction was obtained in the ENP case. The sign of the correction 
is important because a positive sign (as in the case of DL) makes the entropy a concave 
function of the area whereas it is a convex function if the sign is negative or if there is no 
correction (as in the case of ENP). Here we also work with the grand canonical ensemble, 
but we do not set the chemical potential to zero. We employ the counting scheme originally 
proposed by Ghosh and Mitra (GM) ||, but our methods are applicable to the DL scheme 
as well. Because the chemical potential is not set to zero from the start, the entropy ends 
up being a function of the mean area, A, and the mean number of punctures, N, making it 
clear that the positive sign attached to the logarithmic correction is due to a contribution 
coming from fluctuations in N. In the large N limit, both the "inverse temperature" (by 
which we will mean the variable conjugate to the horizon area) and the chemical potential 
(variable conjugate to N) are functions only of the mean area per puncture and so they are 
related to one another as well. At a constant temperature A becomes proportional to N 
and only in this case does one obtain an entropy-area relation whose dominant contribution 
agrees with the BH law and whose sub-dominant contribution is identical to the results of 



TT| for the DL counting scheme. 



We first look at the grand canonical emsemble for isolated horizons, i.e., without the 
projection constraint, in section II. This will set the stage for later calculations and here we 
will show that fluctuations in N lead to an additional logarithmic correction to the entropy, 
coming with a coefficient +1. In Section III we introduce the projection constraint and 
evaluate the partition function by the saddle point method in the large N (thermodynamic) 
limit. We determine the entropy of a black hole horizon as a function of A and N and 
obtain a logarithmic correction to the BH entropy at constant temperature, but it comes 
with a positive coefficient of +1/2, in agreement with the results of |TTJ for the DL counting 
scheme. Thus the effect of the projection constraint is to reduce the number of configurations 
by subtracting | In iV from the expression for the entropy of an isolated horizon. This was 



3 



also determined in |DJ for the canonical ensemble. In section IV we verify the validity of 
the saddle point method by comparing its results with the exact partition function for an 
equispaced area spectrum. We conclude in Section V with a brief summary of our results 
and a short discussion of the outlook for future work. 



II. THE GRAND CANONICAL ENSEMBLE 

We start with a generic isolated horizon and refrain from applying the projection con- 
straint. We want to construct the grand canonical partition function for this system in the 
GM counting scheme, so let us begin with 

00 N ATI 

5(0, a) = E f\~l + l)^e-( 8 ^-^ , (1) 

N=0nj=0 "J 3 ' j 



(as in [ 17 1) where the first term represents the degeneracy of states with N punctures, rij 
of which contribute area ctj = lp\/j(j + 1) with j G {1/2, 1,3/2, . . .} and (2j + 1) is the 
degeneracy associated with spin j. The Barbero-Immirzi parameter is 7 and the variables 
(3 and a are Lagrange multipliers, which can be thought of as an "inverse temperature" 
conjugate to the area and a "chemical potential" conjugate to the number of punctures 
respectively (it is \i = a/f3 that is conventionally called the chemical potential, but we shall 
simply use a). 

Summing over the punctures, 

stf.«) = E(E(«+D«-^^) =r^m- (2) 

where A (a) = e a is the fugacity, and 

^(/3) = ^(2j + l)e- 8 ^. (3) 
3 

The average occupation number of punctures in a state j will be 

1 dlnH _ A(2j + l)e~ 8w ^ 
(%) ~ ~8^~fa~ ~ T^Yz ' [ } 

and the average total number of punctures is given by 

j 

Non-negativity of (N) means that Xz < 1 and, furthermore, the average area of the horizon 
will be given by (henceforth dropping the angular brackets) 

d In E d , . , . A dz N dz . T d In z 

A -~ = TB M1 - Xz) --— z m--TTp = - N l>F- (6) 
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or 



.4 



d In z(a) 



def 



(7) 

47r 7 /2jV da V ; 

where we set a = AirjPlp. Solving this equation would give a = a (a), i.e., the "temperature" 
depends only on the average area per puncture, a. We recognize this as precisely equation 
(9) in fllTf , therefore the relationship between a and a remains intact as we go from the 
canonical to the grand ensemble and from ([|) we find the fugacity, 

N 

On the other hand from (0) and (P), 

lnH = ln(JV+ 1), 
so the Legendre transform of InS, which is the entropy, becomes 

S(A,N) = lnE + (3A-aN = (N + 1) ln(iV + 1) - NlniV + Naa(a) + Nhxz(a) (9) 
and simplifies, in the limit of large N, to 

S(A, N) « In N + N[aa(a) + In z(a)} 



a [a] 



A 



7T7 4/2 



+ A^ln^(a) +lnA^. 



(10) 



This can be compared with equation (12) of It is the general expression for the entropy 
of an isolated horizon of average area A and carrying N punctures. 

First, we notice that allowing fluctuations in the number of punctures has introduced 
an additional logarithmic term in N. The fact that this term is missing in the canonical 
ensemble [17] is an artifact of holding N fixed. Secondly, the entropy in ( ]TQ|) is not the 
BH entropy in general but a complicated function of A and N. However, an important 
consequence of (0) and @ is that both a and A depend exclusively on a in the large N 
limit therefore, for example, a and A remain fixed during any isothermal process involving 
the isolated horizon. Suppose that at some fixed value of the temperature, ao, or of the 
chemical potential, «o, we find that a(a ) = a then 

.4 



N 



47T7^a 



can be used to eliminate iV in ( |10"1) and gives a reduced entropy function 



S(A) w — 

7T7 



cr + 



In z(a ) 
a 



A A 

— + \n — + const., 
4/2 4/2 



(11) 



;i2) 



which is consistent with the BH law, but only so long as the temperature is held fixed. In 
this isothermal condition, agreement with the BH entropy to leading order then determines 
the Barbero-Immirzi parameter in terms of the temperature as well, according to 



1 

7 = ~ 

7T 



lnz(ao) 



(13) 
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An additional criterion is required in order to fully determine all the parameters. For 
example, a possible condition for determining ctq may be taken to be a vanishing chemical 
potential, as in fTI|j , by which the punctures are effectively treated as a photon gas. The 



condition for zero chemical potential, according to (§), is 

z(a ) « 1, (14) 



as obtained in [[L7J, and it determines ao- Solving (|14"D numerically one finds, as in [[J, that 
(T = 0.861, a = 2.921 and 7 = 0.274. 

If the chemical potential is zero it is important to understand how the punctures get 
"thermalized" . In the analogous situation of a photon gas, the photons are thermalized by 
interactions with the walls of the container. It is possible that such a thermalization occurs 
via dynamical processes near the horizon, but it is also possible that such processes cause a 
to be different from zero. 

Note that in the above arguments it has not been necessary to obtain an explicit solution 
for z(a). However, it is possible to obtain a closed form expression for z(a) in a certain 
approximation that will be useful in a subsequent section, when the projection constraint 
is implemented. In this approximation |20| we replace 2j by / so / G {1,2,,...} and aj — > 
« mi + 1). This is an excellent approximation for large values of I but 



not so good for low values of I, so we introduce errors for low spin punctures. Although 
the mathematical steps to estimate the error in this approximation are discussed in |T7[, it 



has been shown that a coherent state representation of the Schwarzschild horizon [21 



or 

a unitary invariant representation of the area operator |22| lead to precisely this kind of 
spectrum. 

For this equispaced spectrum we obtain 

* <0 V) = ^ (15) 



for the partition function. The zero chemical potential condition in (|14|) then gives <7o = 
0.810, ao = 3.318 and 7 = 0.258, values that are not far from the exact ones quoted above. 



III. PROJECTION CONSTRAINT 

With the projection constraint, the canonical partition function for a fixed number of 



punctures was given in [17 



as 

N 



To obtain the grand canonical partition function, it is necessary to introduce a chemical 
potential and sum over N. Therefore, the projection constraint will lead to a partition 
function of the form 

A 1 



n cr, a 



r — I jr e^vw^) sin ^ + 1 ) j e aN (17) 

N ^ _7r ^ \ 1=1 ^ ^ / 
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and performing the sum over N we arrive at 



~[(T,a) = — 

Z7l 



dk 



t 1 - A (a) YaIi z i( a ) 



sinfc(£+l) 
sin k 



(18) 



where we have set Zi(a) = e _<J v /z ^+ 2 ) and A(a) = e Q . 

It does not seem possible to perform the integral exactly, but it may be approximated 
by steepest descent. First we note that that the integrand is an unimodal, symmetric 
distribution centered at k = so long as the denominator of ( |18D is positive. If A = 1 then 
this is the condition that a > <j = 0.861..., which solves equation (fH]). We then extend the 
integration to the entire real line by the change of variables 



k(x) = 2tan~ 1 0/2), 



under which the integrand becomes 

1 



/(*) 



dk(x)/dx 



(19) 



and finally we assume that the distribution is sufficiently well approximated by a normal 
distribution (centred at x = 0) in the limit of large N and to leading order in N, so we 
replace f(x) by 



The variance is given by 



a' 



/(z)«/(0)exp 

'd 2 In/' 
dx 2 



x 



2a 2 



(20) 



/(0) 



and the partition function gets approximated by the area under a Gaussian, which can be 
readily evaluated as 



[a, a 



27r/(0)a 



1 



tt{1 - Xz(a)}{l + Xb(a)} 



(21) 



where 



( ff )=5>(*)(* + l) 
"2 



z=i 

oo 



6(a) =x;^w 



Z=l 



/3 + 2/ 2 + -/ _ 1 

3 3 



(22) 



This gives the average number of punctures as 

<91n5_ X{z(a) -b(a)[l- 2Xz(a)]} 
8a 2[1 + A6(<r)][l - Xz(a)} 

and, once again, the large N limit requires Az(cx) — > V~ . 



(23) 
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If we set 1 — Az(cr) = 5 2 so that 5 2 ~ 0, then to leading order it will be seen that 

5 ^ -^= + 0(N~ 1 ), A = e a « - + ^(iV- 1 ) (24) 



and 



In 5 « - In 2N, 
2 



(25) 



An'yljN da ' 

the last being identical to ([?]). It can be used to determine a — a(a) and z = z(a) as before. 
Putting it all together, one finds that the entropy is now given to leading order in N by 

S(A, N) = - In N + N[aa(a) + In z{a)). (26) 

The only difference between the entropy with the projection constraint and the entropy 
without the projection constraint in (|i~0|) is the coefficient of the logarithmic term, showing 
that the projection constraint has the only effect of subtracting | In N from the entropy. This 



result was obtained in the canonical ensemble as well |L7| where, there being no fluctuations 
permitted in N, a net negative logarithmic correction was found. The positive correction 
in the grand canonical ensemble results from fluctuations in the number of punctures and 
makes the entropy a concave function of the area. 

As before, the entropy is a complicated function of A and N but may be expressed, at 
fixed a = o"o, as a function of the area alone, 



S(A) = — 

7T7 



lnz(oo) 



_i + l ln _i + const . (27) 



For a vanishing chemical potential, the condition z((Tq) w I still holds, according to (p|). 
Thus one also obtains the same value of the Barbero-Immirzi parameter as in the absence 
of the projection constraint. 

As mentioned in the introduction, an interesting approach for the GCE was developed 
in |TTJ for the case of zero chemical potential. The partition functions worked with by the 



authors differ from ours in (|I5D ; in particular, for the DL counting scheme the authors 
evaluate 

s dl (ct) = j_ r ^ — — 

27T J_„ 1 _ 2A(a) COS kl 

with X(a) = 1 and for the ENP scheme, 

Gin ^ h r\)c 

(29) 



, ^E^p/ \ sin k dk 



K J -TV i _ A( a ) E/=i e- CT V^^g±i! 



again with A(a) = 1. The integrand of the first (DL) is also an unimodal, symmetric 
distribution centered at k = and can be treated in a manner identical to ours with the 
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same result ( p6]) in the large N limit. The integrand of the second (ENP) is bimodal, 
symmetric about k — and skew about each mode. The absence of a logarithmic correction 
to the BH law in this case can be qualitatively understood if it can be shown that each mode 
has the effect of subtracting | In N from the entropy. However, our methods are not adapted 
to this distribution and |yj appears to be the most effective approach to this problem at 
present. 



IV. EQUISPACED SPECTRUM 

Given that the distribution function is not in fact Gaussian (if we assume that it is 
Gaussian and compute its kurtosis we should find an excess kurtosis of 3, which would 
contradict our assumption), one may wonder how accurate the saddle point approximation is. 
At present we do not know of any way to formally estimate the error, so we will instead verify 
the leading order discussed above by using the equispaced spectrum mentioned earlier. For 
this spectrum, H can be evaluated exactly. Performing the sum in flHf), with zi(a) = e~°~( l+1 \ 
we find 

1 r dk 
E(a,a) = — - 2cosfc _ e - CT v ( 30 ) 

and the integration may be carried out quite easily; one finds 
where 

A = e 2a + 1, B = 2e a , C = e 2u + \e~ a + 1, D = 2(\ + e a ), (32) 

provided that C/D G (— oo, —1) U (1, oo) for a real E(er, a). Since C and D are both positive 
we require C/D > 1, or 

Xz {0 \a) < 1. (33) 

Finally, using ([32]) and after a little algebra, we obtain the exact partition function for the 
equispaced spectrum, 



Ae CT + V-4(A + e CT ) 2 + (1 + \e.~ a + e 2(T ) 2 

E(a, a) = — ^ =L (34) 

(e Q + e CT ) v/-4(A + e CT ) 2 + (1 + Xe~ a + e 2 -) 2 

and (|33|) will be recognized as the condition that the term under the radical is positive. 

Therefore, using the expression for z^> in (|1"5|), we set, for a positive quantity 5 2 , which 
we do not as yet assume to be small, 

(e ff -l) 2 (2-e- ff ) v ; 

Now the chemical potential may be eliminated in favor or 5 and the partition function 
written in terms of (er, 5). We find 



(1 + 5 2 )(2e* - 1)5 + e 2 *(e- - 1) J 5 ^^ 1)+ ^ 
(a ' 5) = 5[e^ + W° - l)^ 2 ] (36) 
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and taking the derivative with respect to a at constant a, i.e., 

d\nZ(a,a) 1 - 5 2 din Ejo, 5) 

recovers N. 

In the thermodynamic limit we obtain that either e a ~ 1 or 5 ~ 0. For e a ~ 1, we will 
see that a — > oo, which we discard as being unphysical since it implies that N — > for a 
finite A. Therefore we take 5 ~ and consider a small <5 expansion of InS: 



so 



which gives 

and, according to (j35[) 



, n5 = H («--oV^ EH |+ow (38) 



5 « -== + C^iV" 1 ) (40) 



'2JV 



A - { {-^ + ^(iV- 1 ) = -L + O(iV-), (41) 



both of which results agree with (|24|) . Again, if we consider the average area, we find 

O(5 ) (42) 



A dlnE(a,a) 1 - 3e a + 4e 2a e" ff / 2 (l - 3e ff + 4e 2ff ) 



47T7/2 9(7 2(e CT - l)(2e°" - 1)5 2 (e CT - 1) V2e CT - 1 5 

which, to leading order becomes 

def A l-3e< T + 4e 2(T dlnz^ 

a ~ 47r 7 /2iV ~ (e CT -l)(2e CT -l) ~ ' ( ' 



as given in (25) . Solving for a gives 

CT _ 3(1 - a) ± y/( a -l)( a + 7) 
4(2 - a) 



(44) 



Now both solutions are real only so long as a > 1. The solution with the positive sign is 
monotonically increasing from zero to 1/2 with a G [l,oo). This solution is unacceptable 
since e a must be > 1. Thus we conclude that 



u = 3(1 -a)- y/(q-l)(q + 7) 

4(2 -a) 1 j 

and q > 2. As e CT — > oo we see that a — > 2 + and A — > oo. On the other hand, as e a — > 1 + 
then we have a — > oo or N — > and A — > 0. Since both ( [24]) and fl2"ol) are satisfied, the 
entropy will be given by fl26|), with z(cr) replaced by z^°^(cr). 
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V. CONCLUSION AND OUTLOOK 



In this paper we have followed up on work begun in [17|| by providing a detailed analysis 



of the entropy of loop black holes in the grand canonical ensemble (GCE), employing the 
counting scheme of GM. As pointed out in the introduction and in ||12|| , the grand canonical 
ensemble is probably the one that is most relevant to the description of physical processes, 
such as evaporation, which involve changes in the horizon area. 

We observe that in the thermodynamic limit the ratio of the average area to the average 
number of punctures emerges as a fundamental variable that controls both the "temperature" 
(in the area ensemble) and the chemical potential of the black hole. The entropy is in general 
a complicated function of the average number of punctures and the average area, and the 
BH law is recovered only for isothermal processes (equivalently, processes with constant 
chemical potential). 

We have also shown that a positive logarithmic correction to the BH entropy is obtained 
in the GM counting scheme and it is independent of any particular choice of the chemical 
potential. The same can be expected for the DL counting scheme, as we have argued in 
Section III, since it is a unimodal, symmetric distribution centered at k = 0. The source of 
the positive correction was identified as the fluctuations in the number of punctures and the 
projection constraint was shown to decrease the number of configurations by subtracting 
| log A from the entropy for an isolated horizon. The resulting concavity of the entropy as 
a function of the area signals the stability of the black hole in the area ensemble, [Q| , p3[ . 
however such a stability cannot be deemed a "thermal" stability so long as a rigorous relation 
motivated from LQG linking the area of an isolated horizon and some quasi-local energy 
function does not exist. 

An attempt to link the two was made by Ghosh and Perez Jl2| who showed that the 



energy associated with an isolated horizon by a preferred family local observers at a proper 
distance I from it is proportional to its area. The constant of proportionality was determined 
to be the local surface gravity measured by the locally non-rotating, stationary observer and 
Ghosh and Perez were able to work in the microcanonical and canonical energy ensembles. 
Pranzetti |23| extended their construction to the grand canonical ensemble and by matching 



the description of weakly dynamical horizons |25|, £(J with the local statistical description 
was able to obtain a temperature regulating the exchange of energy between the bulk and 
the horizon in terms of the local surface gravity of the horizon. This procedure also involved 
singling out a physical time parameter, equivalently a preferred family of observers, with 
respect to which the evolution of the boundary states was described. 

These works are important steps in the direction of a fully quantum derivation of the 
Hawking radiation from black holes. We feel that rigorously connecting different statistical- 
thermodynamic processes with families of observers via our analysis will give thermo- 
statistical meaning to the choice of vacuum in the semi-classical theory and is the natural 
next step in this analysis. We will report on the results of our findings elsewhere. 
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